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■ ■ Abstract 

Nonlinear cubic Euler-Lagrange equations of motion in the travehng variable are usually derived 

^ ■ from Ginzburg-Landau free energy functionals frequently encountered in several fields of physics. 

t^^ ■ Many authors considered in the past damped versions of such equations with the damping term 

'nI" I added by hand simulating the friction due to the environment. It is known that even in this damped 

t^ ■ case kink solutions can exist. By means of a factorization method, we provide analytic formulas 

for several possible kink solutions of such equations of motion in the undriven and constant field 

driven cases, including the recently introduced Riccati parameter kinks which were not considered 

X 

^ I previously in such a context. The latter parameter controls the delay of the switching stage of 

^ I 

the kinks. The delay is caused by antikink components that are introduced in the structure of the 

solution through this parameter. 
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Nonlinear field excitations occur in a rich variety of collective phenomena. Perhaps 
the most prominent are the fields with cubic nonlinear equations of motion because they 
have effective potential energy with two minima of the type Bip'^ — Aip^ and are commonly 
used as order parameters in the frame of the Ginzburg-Landau theory for the study of 



fl 



'errodistortive domain walls [l|, in the analysis of the phase separation in binary mixtures 
2|, for diamagnetic (Condon) domains |3|, and in more general situations such as domain 



a 



walls in nonequilibrium systems [^. A recently proposed discrete model for the curvature 
modes along protein backbone chains used by Chernodub et al. |5| to explain protein 
folding belongs to the same approach. The dynamics of all these systems can be treated 
variationally as a relaxation process toward one of the stationary states in the potential 
wells. This relaxation is governed by equations of the type dip/dt oc —6J^/6tp where J-" is 
the free energy of these systems. 

In this work, we focus on the solutions of the following equations of motion in the traveling 
coordinate ^ = x — vt and with rescaled coefficients 

^" + pi)' - Bii)^ + Ai^ = (1) 

and 

7/^" + P7/;' - 5i^2 + Ai^ + 7ir/ = (2) 

if a constant external field rj multiplied by its scaled coupling constant 71 is added. Except 
for the friction term, these equations are Euler-Lagrange equations of motions corresponding 
to Ginzburg-Landau functionals. The derivation of such equations can be found for example 

h f 

in Ref. [l| in the context of ferrodistortive domain walls, while recently, Mavromatos [6| 
discussed a well-known counterpart of ([2]) in the case of microtubules where the friction is 
attributed to the so-called ordered water molecules. We make clear that the friction coeffi- 
cient p, although constant, depends on the (constant) velocity of the frame. Both relativistic 
and nonrelativistic dependencies can be encountered in the literature. For example, in the 
case of ferrodistortive materials p oc f (cq — t^^)"'^/^, where cq is the limiting velocity of the 
system, whereas p oc v{D — ?7if ^/2)~^, where D is a diffusion coefficient and m is the inertia 
parameter, in the case of fast, nonoverdamped motion of spin domain walls in Ising ferro- 
magnetics [7|. For equations without the friction term, kink solutions have been well known 
for at least 40 years 8|. In fact, MontroU |8| mentions that Fisher already investigated nu- 
merically solutions of the Fisher equation with a first derivative term. The existence of kink 



solutions in the presence of friction terms has been well settled since the works of Lai 



B, 



Geicke 10|, and Kashcheev ll|. The first goal of this work is to show that for the Eqs. ([T]) 



and ([2]) various kink solutions can be easily obtained by a factorization technique that we 



introduced previously 12|| for equations of the form {Dg = d/ds) 

Dlij + pDs^ + F{4;) = 0, (3) 

where -F(^) is a polynomial function, which in the case ([T]) and ([2]) is a cubic polynomial. 
Equation (^ can be factorized as follows: 

[D,-f2{mDs-fimHs) = o. (4) 



Expanding (^ , one can use the following grouping of terms [12| : 

^'^ - (/i + /2 + ^A Dg^p + f,M = 
and comparing Eq. (|3]) with Eq. ([5]), we get the conditions 



(5) 



/lW/2W = ^, (6) 

^^(^) + ^^ - -^ (^) 

Any factorization like (jl]) of a scalar equation with polynomial nonlinearities of the form 
given in Eq. ([3]) allows us to find a compatible first order nonlinear differential equation, 

[Ds - /i(^)]^ = DstP - /i(V)V' = , (8) 

whose solution provides a particular solution of ([3]). In other words, if by some means we 
are able to find a couple of functions fii'ip) and f2{'ip) such that they factorize Eq. ([3]) in 
the form (jlj), solving Eq. ([8]) allows to get particular solutions of ([3]). The advantage of 
this factorization is that the two unknown functions fi{ip) and f2{'ip) can be found easily by 
factoring the polynomial expression (Q in terms of linear combinations in rational powers of 
ip. This technique is used in the following to find kink solutions in the undriven case and the 
constant field driven case. We also discuss the kinks based on the general Riccati solution, 
which depend on a control parameter of the switching features. The latter kinks, which we 
call Riccati parameter kinks, have not been discussed previously in the Ginzburg-Landau 
framework and drawing the attention to them is another important motivation for this work. 



We first consider the case of zero external field. Montroll showed that Eq. ([T]) has a 
unique bounded (kinklike) solution of the form 



i^MiO = « 



V2: 



a 



(9) 



1 + exp(a^) 

where a = (b — a)/\/2 and the parameters a and b are two of the solutions of the cubic 
equation 

{i/j - a){ilj - b){ilj - d) = ij^ - ij . (10) 

Taking a = 0, b = 1 and d = —1 gives a = 4=. The Montroll kink ipM can be easily derived 
through the factorization procedure just described. Indeed, Eq. ([1]) can be factorized in the 
following two forms (-D^ = ^) 



and 



^ = 



^ = 0. 



(11) 



(12) 



However, factorizations ( TTTl) and (IT2l) are only possible for p± = ±^^^/A^ as obtained from 
Equation f lTT]) is compatible with the Riccati equations 



i)' ± 2^2 (v^^^ - y^itp"^) = 0, 



(13) 



whose (particular) solutions are 



^1,2 



^1 



ki 



,tVAI(o/V2 



On the other hand, the compatible Riccati equations for (IT2l) are 



(14) 



iIj' ± 2-2 (v^A^^ + v^^^) = 



(15) 



with the particular solutions 

^3,4 = - 



^1 ^fA[ 



B^ + QTVM{i-^o)/V2 



fBl + l-eTVMi/V2 ■ 



(16) 



All these solutions are similar to those given by Geicke 10| for this case and one can notice 
that ipi corresponds to the Montroll kink of parameters (0, 1, —1) when Ai = 1 and Bi = 1. 



Moving now to the more complicated case given by Eq. ([2]), let 71 r/ = Aie — Bie^ and 
Lp = il) + e (see also [131]). Then, we get 



^p" + py?' - ^ {Bi'/ - ?,Bie^ + {Wie^ - Ai)) = . 



(17) 



Case I. The factorization of flT7|) can be achieved with /i = aii ^^/Bl^p — ^"^(e)] and 
/2 = (i^^i [VBiif — r_(e)] , where 

3/B7e ± (4Ai - 3Bie^)^ _ S/B^e ± V^T 



r±^ej 



(18) 



From the second factorization condition ([7]), we get 



ai = ±2 2^ 



A-' 



.(-) 



±2" 2 r.fe 



'A, 



(19) 



To have real values of the parameter pj_ , one requires e < AAi/3Bi, and one gets a 



positive-valued friction parameter if e G 



' Ai 2 / Ai 



for the positive front sign of p and 



e e 



"T^a/b ' \/ B ) ^"^^ ^^^ negative front sign of p. We are led to the following Riccati 



equations: 



'B, 



Bi 



^i^Xh^^ T'\/^r+(e)v? = 



(20) 



having as particular solutions 



^t 



1 



2r+(e) 



r-u e 



«! 



where ^0 is a constant of integration. The last step is to go back to the ijj solution, ijjf 



(21) 



^ 



± 



e. Solutions fl2T]) were first found by Geicke through an ansatz method [1 



Case 11. The factorization of (IT7j) can be also achieved with /( = a^/2 and /g = cti^/i. 
Using condition ([7]), one gets 



ai = ±2 2 2 



.(+) 



pL*' 



±2" 2 



r^ie 



'A. 



(22) 



Real values of the parameter pj_ are again obtained for e^ < |^, which implies a positive- 
valued friction parameter if e G {~\/^^^\/ ^ fo^ the positive front sign of p and if 



e G 



-7=-! /^gi, — w -g^ j for the negative front sign of p. 
This factorization implies Riccati equations of the form: 



(^5±2-2(^2^2-2r_(e)(^ = 



(23) 



whose particular solutions are 



vt 



1 



2r_ 



02 



(24) 



/B72 + eT"2(5-&) ' ^^ ^2 

The factorization method allows even more general solutions, the so-called Riccati pa- 
rameter solutions introduced by Reyes and Rosu [ij] which are based on the general Riccati 
solution. Indeed, all Riccati equations in this Brief Report are of constant coefficients, say 
y' — C\y^ — C2y = 0, and if a particular solution yi is known then the general solution 
depending on a free parameter denoted by A can be written as 



y\. 



,Cl,C2 



yi + 



(25) 



A - cih 

where /i = L (2ciyi + C2) dx and I2 = Jc e^^^^'dx. For the nonlinear equations of motion 
discussed here the formulas for the kinks given by Eq. (!25|l are as follows. 

(a) Zero field (kinks having mixtures of rising and decaying exponentials, both of width 

a-i = v/2A4r): 



< 



"M 



T^i 



oc^in-io) 



1 



(26) 



where the minus sign corresponds to the first factorization and the plus sign corresponds to 
the second one. 



^1,; 



^1 



1 + 



1 



\y/A[{l T /Bre"«-«o)) - 1 



(27) 



with the same rule of signs, 
(b) Driving constant field. 
Case I. (kinks having mixtures of rising and decaying exponentials both of width 

ar' = v^/r+(e)): 



and 



^\ 



,,+ _ 2r^(^) 




^A 2 + e-°i(5- 


-?o) 


2r+(6) 


r 


1 


2 + e+"i(^"€o) 


.\ 


/B, 



+ 



1 1 

_l 






v^ 2Ar+(e) [1 + 2e+°i(5-?") - 


V^J 




1 

1 





Bi 2Ar+(e) [1 + 2e-"i«-««)] + 2v^e-"i«-€o)_ 



(28) 



(29) 



These kinks do not have singularities if A ^ fo, ^^^ and A ^ ""'2^''-')' respectively. The 
parametric solutions if)^ are obtained immediately by downshifting the above solutions by 
e. Plots of solutions tljt and ip^ are displayed in Figs. 1 and 2, respectively. 



Case II. (kinks with mixtures of rising and decaying exponentials both of width ag ^ 



The parametric solutions turn out to be 



^t 



2r_ 



2 + e-°2(«-«o) 



1 



^ 2Ar_(e) [1 + 2e"2(5-6j)] _ /g^ 



and 



2r_ 



"Px 



2 + e+°2(?-?o) 



^ 2Ar_ (e) [1 + 2e-°2(«-«o)] + 2v^e-"2«-«o) _ 



(30) 



(31) 



Choosing \^[0,^ 

the forbidden interval for A is 



the (f^ kink does not have any singularities, while in the case of (f^^ 
— ,0). Again, the corresponding tjj parametric solutions 



2r_ 



are obtained by downshifting by e. Plots of the solutions ijj^ and ip^ in this case are displayed 
in Figs. 3 and 4, respectively. 

We now briefly comment on the stability of the Riccati parameter kinks. The stability 
analysis depends on the form of the velocity dependence of the friction-like parameter p. 
If we take this dependence as in the paper of Collins et al. [l|, one can follow step by 
step the stability procedure presented therein. This is because at the first step of the 
stability analysis, that of writing a perturbed kink solution y{C,t) = y{(; A) + 6y{(,t) one 
notices that the Riccati kink y{(] X) is a solution of the same equation as the common kink 
y{(; A = oo) and therefore the linearization leads to the same eigenvalue problem. For a 
velocity dependence of p corresponding to a driven, damped, nonlinear Klein-Gordon type 
equation, the stability analysis is somewhat more complicated but has been sketched in the 
important paper of Biittiker and Thomas [li 



In summary, using the factorization method introduced in 12|, we have obtained the 
analytic forms of various kink solutions of the nonlinear cubic Euler- Lagrange equations of 
the damped type in the traveling variable. The results presented here can be directly applied 



to the Condon domains if we make t 
parameters given by Gordon et al. 



le following identification of our parameters with the 
sl: p = v/Kr,Ai = A/K,Bi = B/K,-fi = a/k; in 



the case of Collins et al. [l| one should take Ai = Bi = I. Examining the formulas for 
the Riccati parameter kinks, one can easily infer that the parameter A occurs as a control 
parameter of the initiation of the switching stage [ij]. Indeed, A is associated with the 
exponentials of opposite exponent in the denominators with respect to the exponential of 
the particular Riccati solution and this is what generates the delay. Interestingly, we notice 



that although the switching is more delayed when A increases, this is so only at relatively 
low values of the parameter, while at higher values of A the delay saturates. Since switching 
is related to microscopic restructuring of the kink (mesoscopic domain), one may think 
that the A parameter can characterize the dependence of the switching delay on the rate 
at which an applied field is ramped up or down. Finally, all the kinks discussed here occur 



E 



in conditions of environmental friction, which is not easy to define microscopically [6|. If 
the frictional effects are considered as first derivative terms in cubic nonlinear equations of 
motion as done here, then the kinks discussed in this work occur only for very particular 
values of the friction coefficient given by p±, p^T , and pj_ . 
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FIG. 1: Case I. Plot of ^+(^) for A = 0.125 (solid line), 0.2 (long-dashed line), 0.5 (dashed line) 
and 10 (dotted line). Ai = 3, Bi = 0.7, p = 0.90326 (e = 2.2772), and ^o = 0. 
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FIG. 2: Case I. Plot of V'a (0 for A = 0.01 (solid line), 0.1 (long-dashed line), 0.5 (dashed line) and 
10 (dotted line). Ai = 3, Bi = 0.7, p = 2.39335 (e = 1.0351), and ^o = 0. 
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FIG. 3: Case II. Plot of Va (?) for A = 0.77 (solid line), 0.9 (long-dashed line), 2 (dashed hne) and 
10 (dotted line). Ai = 0.7, Bi = 3, p= 1.51635 (e = 0.5313), and ^o = 0. 




FIG. 4: Case II. Plot of tp^ (^) for A = 0.53 (solid line), 0.6 (long-dashed line), 1 (dashed line) and 
10 (dotted line). Ai = 0.7, Bi = 3, p = 0.435766 (e = -0.5313), and Co = 0. 
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